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A hypothetical violation of Lorentz invariance in the electrons' equation of motion (expressed 
within the Lorentz-violating extension of the standard model) leads to a change of the geometry 
of crystals and thus shifts the resonance frequency of an electromagnetic cavity. This allows ex- 
perimental tests of Lorentz invariance of the electron sector of the standard model. The material 
dependence of the effect allows to separate it from an additional shift caused by Lorentz violation in 
electrodynamics, and to place independent limits on both effects. From present experiments, upper 
limits on Lorentz violation in the electrons' kinetic energy term are deduced. 
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INTRODUCTION 



Special relativity and the principle of Lorentz invari- 
ance describe how the concepts of space and time have 
to be applied when describing physical phenomena in flat 
space-time. Improving the accuracy of the experimental 
verification of these fundamental concepts is of great in- 
terest, also because a violation of Lorentz invariance is be 
a feature of many current models for a quantum theory 
of gravity, e.g., string theory 0,|j|, loop gravity 0,01, 
and non-commutative geometry Such a violation of 
Lorentz invariance is described in the general standard 
model extension (SME) developed by Colladay and Kost- 
elecky 0. According to it, Lorentz violating terms might 
enter the equations of motion of bosons and fermions. At 
first sight, the quantum gravity induced corrections and 
effects are of order E/Eqq ~ 1CU 28 where E is the en- 
ergy scale of the experiment (which in ordinary optical 
experiments is of the order 1 eV) and Eqq ~ -Epianck is 
of the order of the Planck energy. Therefore these effects 
seem to be far from being observable in laboratory ex- 
periments. However, as it occurs e.g. in scenarios lead- 
ing to a modification of the Newton potential at small 
distances, some mechanism may apply which effectively 
leads to much larger effects in the laboratory. It is thus 
interesting to find experimental configurations in the lab- 
oratory that can place strong upper limits on as many of 
these terms as possible. 

Experiments on Lorentz symmetry that study light 
propagation have a long and fascinating history, starting 
from the original interferometer experiments of Michel- 
son in Potsdam and Michelson and Morley J8| in 
Cleveland. Modern versions of these experiments |tA lid 
IH mill 111 111 replace the interferometer by a mca- 
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surement of the resonance frequency v = mc/(2L) of an 
electromagnetic (Fabry-Perot) cavity, that is given by the 
velocity of light c along the cavity axis, the cavity length 
L, and a constant mode number m. Lorentz violation 
causing a shift of c or L connected to a rotation or boost 
of the cavity frame of reference can thus be detected 
through the corresponding shift of the resonance fre- 
quency. From such experiments, upper limits on a tensor 
(&f) k\^i> have been found, that encodes Lorentz violation 
in the photonic sector of the SME [11 mm 0013 • 
These experiments are mainly based on the shift of c con- 
nected to non-zero values of {kp) K \^ v , as an additional 
change of L caused by {kp)^^ and a corresponding ori- 
entation dependent modification of the Coulomb poten- 
tial is negligible for most cavity materials [lflj . 

In this work, we treat the effect of Lorentz violation 
within the fermionic sector of the standard model exten- 
sion in cavity experiments. A modified kinetic energy 
term entering the non-relativistic Schrodinger hamilto- 
nian of the free electron (p 2 + 2E'- k PjPk) / '(2m) leads to 
a change of the geometry of crystals, and thus a change 
8 e -v of the resonance frequency of a cavity made from 
this crystal. Here, pj is the 3-momentum, m the electron 
mass, and E'- k = —Cjk — coo^jfe a dimensionless 3x3 ma- 
trix given by a tensor of the SME. Thus, the total 
shift of the resonance frequency Sv = 5 e -v res + ^EM^res, 
where ^em^ denotes the shift due to Lorentz violation in 
the electromagnetic sector. Since 5 e -v depends on the 
cavity material, it can be be distinguished from Lorentz 
violation in electrodynamics by comparing cavities made 
from different materials. Experiments using suitable con- 
figurations of cavities can place separate upper limits 
on the components of e^„ and (/cf)kA/w Using data 
available from past experiments, we deduce approximate 
bounds on some combinations of components of at 
the 10~ 14 level. From future cavity experiments, Earth 
and space-based, that are projected as tests of Lorentz vi- 
olation in electrodynamics |20l l2ll l22j , bounds at a level 
down to 10 -18 can be expected. 
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In Sec. [HI we present the non-relativistic hamiltonian 
of the free electron within the SME. Since violations of 
Lorentz invariance are certainly small, it is sufficient to 
work to first order in the Lorentz-violating modifications 
throughout. The crystal adjusts its geometry such as to 
minimize the total energy (Sh) + C/ C iast, where (Sh) is the 
expectation value of the Lorentz-violating part of the 
hamiltonian and t/ c iast is the elastic energy associated 
with distortion of the crystal. It is calculated in Sec. 
IIIII The resulting geometry change is calculated in Sec. 
IIVI A fairly detailed model for the crystal allows us to 
obtain specific results for practical materials, including 
quartz and sapphire. In Sec. [V] we discuss experimental 
configurations and obtain bounds on Lorentz violation 
in the electrons' equation of motion from present exper- 
iments. In appendix ^ we discuss the hypothetical case 
of a cavity made from a spin-polarized solid, which al- 
lows to place experimental limits on an additional, spin 
dependent term from the SME, at least in principle. In 
appendix [Bj we summarize some conventions made in 
elasticity theory that are needed for our calculations, and 
in appendix [CI we give in detail the Fourier components 
of the signal for Lorentz violation in laboratory experi- 
ments on Earth. 



m e and mP are the electron and the proton mass, 7^,75 
and <j^ v are the conventional Dirac matrices, and D v is 
the covariant derivative. In this work, we deal mostly 
with electrons and add a superscript to denote param- 
eters for particles other than the electron, e.g., c M „ is 
a parameter for the electron and d^ v the corresponding 
parameter for the proton. The tensors entering M have 
the dimension mass, the others are dimensionless. H^ u 
is antisymmetric; gx^v is antisymmetric in its first two 
indices. c M „ and are traceless. Gauge invariance and 
renormalizability excludes e v ,f„, and gx^, so these are 
either zero or suppressed relative to the others |2^| . 

Lorentz violation for the photons is encoded in the ten- 
sors {UafY and {k F ) K x^- The four degrees of freedom 
contained in (£^i?) K are constrained strongly in measure- 
ments of cosmological birefringence [TtI Il8| and are ne- 
glected in what follows. 10 of the 19 degrees of freedom 
of ( kF) K ,Xnv are constrained by astrophysical observations 
[17L Hal, the other 9 can be measured in cavity experi- 

mentell! H H [H El . 



B. Modified non-relativistic hamiltonian 



Free electron 



II. STANDARD MODEL EXTENSION 
A. Model 

The SME starts from a Lagrangian formulation of 
the standard model, adding all possible observer Lorentz 
scalars that can be formed from the known particles and 
Lorentz tensors. Taken from the full SME that contains 
all known particles, the Lagrangian involving the Dirac 
fields ijj e of the electron and i/j p of the proton and the 
electromagnetic field F^ v can be written as (in this sec- 
tion, we use units with h = c = 1; the greek indices take 
the values 0,1,2,3) HE1 



(1) 



-\(k F ) K x^F KX F^ + \{k AF Te KX ^A x F^ , 

where h.c. denotes the hermitian conjugate of the previ- 
ous terms, and A x is the vector potential. The symbols 
r®' p and M e,p are given by 

T,y = 7„ + + d M v757 p + e v + ifu'fs + ^gx^^ Xfl , 



(2) 



where the superscripts e and p are to be added to the 
symbols a M , 6 M , c^, d^, e M , / M , g Xll v, and that repre- 
sent tensors encoding Lorentz violation for the fermions. 



The non-relativistic Schrodinger hamiltonian h = h + 
Sh of a single free electron within the SME derived from 
this Lagrangian (using Foldy-Wouthuysen methods, |24() 
is the sum of the usual free-particle Hamiltonian h and a 
Lorentz-violating term [23, UM 
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with the components of the 3-momentum pj and of the 
Pauli matrices a 3 . The latin indices take the values 
1,2,3. (We denote both 3-vectors such as Xj and recip- 
rocal 3- vectors such as pj by subscript^ A Hamiltonian 
of this form has also been derived in |2f|. The constant 
term mA' has no physical consequences and is included 
for completeness. The term proportional to Cj can be 
eliminated by choosing coordinates such that the systems 
centre of mass is at rest [25|. The dimensionless coeffi- 
cients A' , B'j,C'j, D'j k , E'j k and F'- kl can be expressed in 
terms of the quantities entering the Lagrangian [2^, |^[ : 
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2. Interaction terms 

In addition to 6h, the Hamiltonian arising from 
the Lagrangian Eq. |QJ also involves modified 
interaction terms proportional to combinations of 
a^b^c^jd^^^f^gx^, and H^. For the non- 
relativistic electrons in solids, however, these are sup- 
pressed by a factor given by a, the fine-structure con- 
stant, relative to the modifications of the free-particle 
Hamiltoninian, Eq. |2^|. This is basically because 
the typical energy scale for such electrons is the Rydberg 
energy ma 2 /2. We can therefore neglect the modified 
interaction terms. 



C. Coordinate and field redefinitions 

Some of the parameters contained in either the photon, 
electron, or proton sectors of the Lagrangian Eq. can 
be absorbed into the other sectors by coordinate and field 
redefinitions without loss of generality. Thus, not all of 
the coefficients in the Lagrangian have separate physical 
meanings. Loosely speaking, in experiments where one 
compares the sectors against each other only differential 
effects arc meaningful. 

For example, in a hypothetical world containing only 
photons and electrons, the nine components of {kF) K \iiv 
not constrained by astrophysical experiments could be 
absorbed into the nine symmetric components of c^ v 
H E3, EE EH- By definition, either the photon or 
the electron sector could be taken as conventional with 
all the Lorentz violation in the other sector. For ex- 
ample, for tests of Lorentz violation for the photon 
[l^llaliall7lll8f . one implicitly assumes a conventional 
electron sector. 

The presence of protons (and neutrons) in the solid 
changes this picture. We can still assume that one of 
the sectors is conventional, but then in general the other 
sectors are Lorentz-violating. Choosing a conventional 
proton sector allows us to disregard the proton terms. It 
also fixes the definition of coordinates and fields so that 
the components of cannot be absorbed into {kF) K \iiv 
in general, i.e. they acquire separate physical meanings. 

To illustrate this, it suffices to consider an extension 
of the toy version of the SME introduced in ^3 • Its La- 
grangian describes electrons and protons as scalar fields 
(jf and <jjP, neglecting spin effects: 



+ (L> AI <j!> p ) t L>^ p - (m p ) 2 (0 p ) t p 
- \f^F^ - \(k F ).x^F^F^ . 



(10) 



It has a conventional proton sector and non-conventional 
electron and photon sectors. Lorentz violation for the 
electron is given by the coefficient k^u- As usual, the co- 
variant derivative is given by D tl (j> e ' p = (d tl +iq e ' p A M )0 e ' p , 
where q e,p is the particle's electric charge. If one iden- 
tifies k^v = c^ v , this Lagrangian leads to the modified 
Hamiltonian Fq. J3J) if only the c M „ are nonzero. 

For simplicity, consider the special case of only one 
nonzero component fc o = k 2 — I , where k deviates slightly 
from unity. The lagrangian takes the form 0] 



{D^D v <j> e + (k 2 - I)|A^ e | 2 ~ m 2 ^ 6 )V 



+ (L>, I </> p ) t L"> p - (m p ) 2 (0 p ) V 
-\f^F^ - ~(k F ) K x^F KX F^ . 



(11) 



By coordinate transformations t — > tk, x — > x, the field 
redefinition Aq — > Aq, Aj — ► Aik, and rescaling the elec- 
tric charge q — ► q/k, one obtains the Lagrangian 

+ (k- 2 - l)\D r\ 2 - (mP) 2 (r)U p (12) 



1 



-B 



where B is the magnetic field. Thus, the Lorentz viola- 
tion in the electron sector has been moved to the proton 
and photon sectors. (If (kp)^^ = 0, the parameter 1/fc 
can be interpreted as a modified velocity of light 
However, it is in general not possible to eliminate Lorentz 
violation in more than one sector at the same time. 

Cavity experiments compare the velocity of a light 
wave to a length defined by a crystal. In the light of Eq. 
(|f f \ . the Lorentz violation for the electron acts via the 
term (k 2 — I )\DQ(f> e \ 2 . With a time-independent Coulomb 
potential Aq = const., this contributes a term 



2im(k 2 - 1)[4> i0 + qSR(A )4>} 



(13) 



(V„u + k MU )(D ll <j) e y D v <f> e - m 2 m^ e 



to the equation of motion for <i> e = e~ tmx °(jf in the 
non-relativistic limit (obtained, in the usual way, by the 
Euler-Lagrange equations and setting to zero terms of 
order |^4o| 2 and m°). 3? denotes the real part. The sec- 
ond term modifies the coupling of the electron to the 
Coulomb potential, causing a geometry change of the 
crystal. Thus, a combination of k 2 — 1 and the modi- 
fied velocity of light given by Qsf)kXhv is measured in 
the experiment. 

In the alternative description by Eq. 1)120. the same 
Lorentz violation acts via a term analogous to Eq. 113|) 
in the equation for the proton, i.e., a rescaled coupling of 
the proton to the electric field, and a modified velocity 
of light given by (fcp) K A^ and k 2 — 1. Physically, both 
pictures are equivalent. 
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Here, we considered only a single parameter analogous 
to coo, that causes a scaling of the solid that is rotation- 
invariant. Thus, it cannot be measured in usual cavity 
experiments, that search for a modulation of the effect 
connected to a rotation of the cavity in space. How- 
ever, the tensors c M1 , emerging from our special case via 
the three Lorentz boosts can — i.e., at least three out 
of nine degrees of freedom contained in the symmetric 
part of c^ v . It is not impossible that, by coordinate and 
field redefinitions, some of the other parameters can be 
absorbed into quantities that have no measurable effect. 
However, as we have shown, at least 12 out of 18 pa- 
rameters from the photon and electron sector (restrict- 
ing the electron sector to those parameters that are not 
constrained by cosmological experiments) would be sep- 
arately measurable, that can e.g. be chosen as three c^ v 
and nine {k F ) K \^. 

In what follows, we adopt a conventional proton sector, 
with all the Lorentz violation in the electron and photon 
sectors. One could possibly extract the measurable quan- 
tities from c^ v and consider only those in what follows; 
however, thereby one would single out a preferred frame 
in which the measurable quantities are defined, and loose 
the covariance under observer Lorentz transformations 
which otherwise holds in the SME. Therefore, we choose 
not to do so and treat all elements of as independent. 



D. Previous experimental limits on electron 
parameters 

It is convenient to express limits on the coefficients 
within a sun-centered celestial equatorial reference frame 
as defined in . The components of quantities given in 
that frame are denoted by capital indices. Limits for 
many particles, including muons . protons, and neutrons, 
have been studied, see and references therein. 

For the electron, the limits given below have been found. 
However, to our knowledge there are no experimental 
limits on E'- k and on many components of F'^ kl for the 
electron. 

From clock comparison experiments |23| . a limit on 
Bj < 1CT 24 {mBj is denoted bj in|2|]) is obtained. Fur- 
thermore, for the linear combinations 



dj = m(d 0J + d JQ ) 



md 



1 



JO 



9d,j 



mejKL{gKOL + t;9klo) - bj 



(14) 



dj/m < 10~ 19 and gD,j/m % 10 -19 . These are order-of- 
magnitude limits, since some assumptions are needed to 
extract them from the measurements |23j . 

An experiment using spin polarized solids yielded 



\B> 



(2.7 ± 1.6) x 10 



-25 



in a similar ex- 



and 



periment {{B' x f + (B' Y ) 2 ) 1/2 < 6.0 x 10" 

\B' Z \ < 1.4 x 10" 25 have been found. 



2(1 



Hydrogen spectroscopy could prospectively limit linear 



the superscripts p denotes parameters for the proton) to 
about 10 -27 GeV [3l]]. Comparing the frequencies of hy- 
drogen masers, 32] find \BjmP + B'jm\ < 2x 10~ 28 GeV 
(m and m p are the electron and proton mass, respec- 
tively) . 

The potential for further tests of Lorentz invariance 
in space is discussed in 33]; for the electron, limits on 
several parameter combinations are expected using 1 Cs 
and 87 Rb clocks. However, these tests allow no limits on 
the components of E'j k . 

III. INFLUENCE OF Sh ON SOLIDS 

To first order in the changes, the influence of Lorentz 
violation in the electron's equation of motion on the prop- 
erties of a crystal is induced by the expectation values 
of the Lorentz-violating contribution to the hamiltonian, 
that is calculated in this section. In Sec. 1111 Al we present 
our ansatz for the electron wave function; the expectation 
value (6h) is then calculated in Sec. IIII Bl 

We denote by (x a )i, {p a )i, and (a a ) 1 the spatial, 3- 
momentum, and Pauli matrices for the ath particle. The 
non-relativistic single-particle hamiltonian for the a th 
particle is denoted h a = h a + Sh a . The hamiltonian of 
the solid 



tall 



a^b 



15) 



is the sum of h a + 5h a over all particles, plus the sum of 
the interaction terms h a ,b over all pairs, and over <5/i aj &, 
a possible Lorentz-violating correction to it. (The factor 
\ corrects for the double-counting of pairs.) The Lorentz 
violating terms are contained in Sh a and Sh a ,b- To first 
order in the changes, the resulting modifications of the 
properties of the solid are the sum of the modifications 
arising form the individual terms. 

The interactions in a solid are electromagnetic. The 
geometry change of crystals as a consequence of the mod- 
ification to the interaction term form the photonic sector 
of the SME [rH Irsf has been treated for ionic crystals in 
[l^] . We will not consider this term further here. In this 
work, we deal with the modifications due to the Lorentz 
violation in the electrons' equation of motion, ^ a Sh a . 



combinations of B'j, B f, djo, d p J0 ,HjK, and H P JK 



A. Wave function ansatz for the solid 

According to the Bloch theorem (|34(, pp. 133-141), 
the single-electron wave function tp a for the a-th electron 
(a = 1 , . . . , N) of a solid can be written as the prod- 
uct of a plane wave exp{j-q a x a } (where q a is the quasi- 
momentum of the a th electron) and a function u^{r) with 
the period of the lattice. u^ a (r) depends on q a , and thus 
on the electron number a. To make a Fourier expansion 
of u^ a (r) , we note that, if kji is the 3x3 matrix containing 
the primitive reciprocal lattice vectors ki, any reciprocal 
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lattice vector can be expressed as a linear combination 
riikji with some coefficients «j e Z (HH, pp. 86-87). 
Therefore, 

"0a = —7= exp I jVaXa \ ^{Cn)a exp{-mj ifcy (x a )i} , 

n 

where V is the volume of the (macroscopic) solid con- 
sidered. The (cft) a are the Fourier coefficients of u^ a (r); 
they depend only on rijkij — (q a )i, i.e., they can be ex- 



pressed as c ri 



ffcij-(9a)i 



The r? summation is carried out 



over 1? . The (c n ) a satisfy 



Ei( c »)-i 2 = 1 



(17) 



because of normalization ((VOaKVOa) = 1- If we assume 
that the crystal has inversion symmetry, the origin of the 
coordinate system may be chosen such that f |34j. p. 137) 



(C-n)a — ( c rl)a 



(18) 



The star denotes complex conjugation. 

The normalized antisymmetric TV-electron state $! N 
(N is the total number of electrons) can be constructed 
from the N x N matrix 



^i(fi) ... 4>i(x N ) 



ip N {xi) ■ ■ . iPn(x n ) 



(19) 



as the Slater determinant [35[ ^f N — -^det(tp) [det de- 
notes the determinant of a square matrix] . 

B. Calculation of matrix elements 

1. Specifications 

For a bound system in its rest frame, the expectation 
value of the particle momenta vanishes 



£<(p«)i> = o. 



(20) 



vanishes. The case of spin polarization will be considered 
in the appendix. Furthermore, we assume a vanishing 
sum of the helicities of the electrons, 



and 



£<ba),K) fc ) = o 



£<(p«MPa)i(ff„) k >=0. 



(22) 



(23) 



Although situations could be imagined which violate the 
last two relations in spite of S l = 0, this can be considered 
highly unrealistic. Therefore, the expectation values of 



the terms proportional to Cj, D 



and Fj kl from the 



hamiltonian, Eq. (JSJ), vanish. Disregarding the constant 
term proportional to A, 



{(Sh)) 



m £• — * 



2. Calculation of 5h 



We now calculate the matrix element ((j>i)i(pi)j) for 
the first electron. Since it turns out to be independent 
from the electrons number (a consequence of the anti- 
symmetry of the iV-electron state *f? N ), the sum of the 
matrix elements for all electrons can then be obtained by 
multiplying with the total number of elec- 

trons N. We have 



{(piUpih) = I * 

Iv 



N* 



-h 2 d 2 



^ N d 3 x 1 . . . d*X N . 



d(xi)i(a;i) 



( 24 ) 

The integrations are carried out over the volume V of 
the solid. \E ,Ar and *ff N * are given by Slater determinants. 
Evaluation of the matrix element starts by an expansion 
of these determinants with respect to the first column, 



N 



det(V>)= £(-1)^(1^°), 



(25) 



We also assume no spin polarization, i.e., the sum of the 
spin expectation values 



(21) 



where ) denotes the (N — 1) x (N— 1) minor matrix 

obtained from the N x N matrix (ip) by deleting the first 
column and the a-th row. The derivatives can then be 
carried out: 



1 N 

((j>l)i(pi)j) = ^yj £ (~ 1 ) a+b E( c " l ^( c «)b [(9b)i(Qb)j + h 2 nikun k k jk - hn k hk{qb)j ~ finikji(q b ) 



(26) 



„b=l 
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V 



exp< i(xi)i 



(rrii - ni)ku + z [(©>)* - (la)i] 



d 3 xi / det(|iV'* Q )det(|i^ )d 3 x 2 ■ ■ ■ d 3 x N . 



Note that only the first integral in Eq. (|26() contains x\. 
Using the abbreviation 

Ki = (m l - ni)k u + -A{<lb)l - [q a )l] , (27) 
the d 3 a;i-integration in Eq. I|26|l can be expressed as 

Since the quasi- momenta (q a )i are within the first 
Brillouin zone (0, p. 89), |(g a )/| < (h/2)\n j ki J \ for 
any rij G {Z\0}. Thus, ki = only if n, = wii and 
(lb)l = (la) i- We may assume that this holds only for 
a = b. That allows to carry out the rn and the b summa- 
tions. We now use 

^2\(cn)a\ 2 ni = 0, (29) 

n 

which follows from Eq. Ijl8(l and eliminates the terms 
linear in n from the first line of Eq. (|26[1 36] . We define 

(£,a)ik ■= ^ \ ( c n)a\ 2 nin k (30) 

n 

and 

M7 : = JJ^2(<la)t(qa)j , (31) 

a 

the average of the quasi-momentum product qiQj over all 
electrons. We also define the average 

~: jr ^ iL h ■ (32) 

a 

Together with Eq. (|T7J> . this yields 

((Pi)i(Pi)j) = J7j [WqJ+ff&khikjk] 

x / d 3 x 2 ■ ■ ■d 3 x N det(\ 1 iJ* a )det(\i'$ b ) (33) 
Jv 

To prove that ^\J/ Ar_1 |\p Ar_1 ^ = 1 we expand the remain- 
ing determinants of (|i'0* ) in terms of the column which 
is now the first one. The d 3 x 2 integration can then be 
carried out. The procedure is repeated, until the d 3 x^ 
integration is done. Each step reduces the dimension 
of the Slater determinant by one and produces a factor 
equal to the number of electrons still involved. Taking 
all these factors together cancels the normalization factor 
1/(N — 1)!. Thus, we obtain the desired result 

(ipi)i(pi)j) =WH +f^ilkkakjk- (34) 

Since the right hand side of this result contains no refer- 
ence to the electron's number, it holds for all o = 1, . . . , TV 
electrons. 



3. Estimating 

The properties of the wave function enter the momen- 
tum expectation value via and qjq]. 

£lk is an ensemble average over all electrons: 

£' fc = Z~^2 ( TV K c ") a l 2ninfe ) =: \cn\ 2 nin k , (35) 
H \ a In 

obtained by substituting Eq. (j30J into Eq. 1(32]). A de- 
tailed evaluation of would start from material specific 
Fourier coefficients eg obtained experimentally or the- 
oretically (see, e.g, (32|). Since detailed wave- function 
calculations for realistic materials are notoriously diffi- 
cult and would have to be performed for each individual 
material, it is interesting to use a relatively simple model 
for the Fourier coefficients. Such a model might already 
be quite accurate, since only the average of the absolute 
square \ca\ 2 is required, rather than the detailed (cn) a 
for the individual electrons. A (possibly complicated) 
dependency of the (c^) a on ft can be hoped to smooth 
out in the averaging. The model must, however, be in 
accordance with the requirement that the wave function 
has the rotational symmetry of the lattice, i.e, Cpn = eg 
if p is any operator of the rotation group of the crystal 

(S3, pp. 469). 

For such a simple model, we assume that |c^| = |cijji|, 
i.e., the average of the absolute squares of the average of 
the Fourier coefficients depends only on \n\. It follows 
that 

ilk = 7mat<5ifc (36) 

with some material dependent constant 7 ma t. For deter- 
mining 7 m at, we note that the average kinetic energy of 
an electron (T) = {PiPj) %■ From Eqs. 
and (EH, 

h 2 

(T) = jmnt-^—SikkiikjkSij . (37) 

If we neglect for the moment the energy of the chemical 
bonding, this should correspond to the average of the ki- 
netic energies of the atoms' valence electrons, which can 
be estimated using the Bohr model. The kinetic energy 
of an electron in the Bohr model TboIu- = EjiZ/n 2 . Here, 
En ~ 13.6 eV is the Rydberg energy, Z the charge num- 
ber of the atom core, and n the principal quantum num- 
ber. Thus, the kinetic energy from the Bohr model, aver- 
aged over the electrons within the atoms of the molecule, 

T^Z = — — J2 ^-^T 1 ( 38 ) 
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where N e ^ m gives the number of valence electrons per 
molecule. The index k enumerates the atoms of the 
molecule and Vk,rik, and Zk are the valence, principal 
quantum number, and charge number, respectively, of 
the atom k. Note that Zk — Ufe, since in an atom with 
v valence electrons, the charge number of the atom cores 
and the inner shell electrons is Z = v. For example, in 
quartz, SiC>2, there is one Si atom with v = 4, n = 3 and 
two O atoms with v — 2 and n — 2. Thus, for example, 

(T B ohr) quartz = -jp + = sq Er ■ ( 39 ) 

Comparing this to Eq. (|37|) , we obtain the material spe- 
cific values 7 ma t given in Tab. [I] 

The model can be refined by taking into account the 
energy of the chemical bonding, which leads to an in- 
crease of the actual kinetic energy of the electrons. The 
so called enthalpy of formation A/iJ° gives the enthalpy 
for the formation of the crystal from the elements in their 
usual state at standard conditions (room temperature 
and pressure), e.g., solid or diatomic (O2, for example). 
The Bohr model, however, predicts the energy of the un- 
bound atoms. That means, the change T— TBohr between 
the sum of the kinetic energy of the valence electrons of 
the free atoms TBohr and the sum of their kinetic en- 
ergy in the molecule T is the difference of A/i7° and the 
applicable enthalpies of sublimation A SUD i-£f or dissoci- 
ation Adiss-ff of the elements. For sapphire, AI2O3, for 
example, 

A f H°(Al 2 3 ) = 16.8 eV 
-A subl ff°(2Al) = -2x3.0eV 
-A diss ff (3O) = -3x2.5eV 

T-Tsohr = 3.3 eV, 

or about 5% of TBohr = 68 eV. This would lead to a 5% 
increase of the factor 7, na t. This indicates that the energy 
of the chemical bonding is, for our purposes, negligible. 

The ultimate refinement of the model would be the 
insertion of material specific Fourier coefficients into 
Eq. H35fl . The precision of the model would then ap- 
proach the limitations of the Bloch ansatz for the wave 
function itself, which is based on a mean field model for 
the electron-electron interactions. Such a detailed model 
is, however, beyond the scope of the present work. 



4- Result 



The quasi- momenta \q a \ are restricted to the first Bril- 
louin zone, \q\ < Most electrons, however, will 

have a quasi-momentum lower than this maximum value, 
so that the average q^qj is a relatively small contribution 
to (8h). If we neglect it, we obtain the final result 

h 2 

(Sh) = N^t—E^Sikkukjk . (41) 

fiikkukjk is symmetric in the indices i and j. This result 
will be used in the next section to compute the geometry 
change of the crystal caused by Lorentz violation in the 
electrons' equation of motion. 

IV. CHANGE OF CRYSTAL GEOMETRY 

The direct lattice vectors l a contained in the matrix 
Ua determine the structure of the lattice without Lorentz 
violation. Lorentz violation will cause a change h a of the 
crystal geometry, i.e., the lattice vectors are now given by 
ha + ha- To calculate it, we adjust h a such as to minimize 
the total energy of the lattice 

U = U (ka + ha) + (5h(lia + ha)) ■ (42) 

The first term is the conventional total energy of the 
lattice without Lorentz violation. It can be expressed as 

Uo{ha + ha) = U (ha) + C4last(L) + U c (43) 

where U c is a constant and C/ C i ast is the elastic energy 
connected to a distortion of the lattice. If h a = 0, the 
elastic energy J7 i as t = 0. The total energy is thus given 
by 

U = Uo(ha) + U clast (h a ) + (5h{l la + l la )) . (44) 

The correction ha is found by setting to zero the deriva- 
tive: 

r)TI BIIi , d {Sh(ha + ha)) 
dljb dljb dljb 

To explicitly calculate Z^, we have to express the contri- 
butions to U in terms of Ijt,. We will do so for (Sh) in 
Sec. HVAl and in Sec. ITVBl for C/ clast . In Sec. WU the 
total energy per unit cell thus obtained is minimized and 
the geometry change as expressed by a strain tensor e,j 
is calculated. 



The expectation value of the Lorentz- violating correc- 
tion to the single particle hamiltonian ^ Q 8h a can be 
obtained from Eq. (|34|l by multiplying with the number 
of electrons. Inserting ^ as obtained in the previous 
section into Eq. (|34|l . we obtain 

(Sh) = —E' i:j (q^qj + -y ma .th 2 Sikkukjk) ■ (40) 



A. Dependence of Sh on U a 

The change of the hamiltonian's expectation value 
(Sh), Eq. ijlTJl depends on the geometry via the recip- 
rocal lattice vectors kij, for which we have the relation 
(0, P- 87), 

hjk ik = 2irS jk , (46) 
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and therefore hjk n j = 2nS in . If we substitute kj + 1^ 
and kjk + kjk, with 1^ <C lij and kjk <C kjk, we obtain 



(kj + hj)(kik + k lk ) = 2nSj k 



(47) 



with some coefficients r\ a g N. Multiplying this equation 
by fcbc and using Eq. (|46[) . we obtain 



?7c — 7: kb c Xb ■ 

Zir 



(56) 



hjkik + hjkik + hjkit + hjkik — zirSjk . (48) 

The first term on the l.h.s. cancels with the r.h.s due Eq. 
(|46|l : we neglect the second order term on the l.h.s., and 
obtain 

kjhk + hjkik = . (49) 
Multiplying with k n j and using Eq. 146(1 again, 



If l{, a is shifted to l ba + lb a , the lattice point originally at 
Xb will be shifted to Xb + Ub, where 



u b — Va^ba — „ ^da^d^ba • 
Z7T 



Therefore, we have 



^-—k I 
dx ~ 2ir ca 



(57) 



(58) 



Z7T 



(50) 



We can now substitute kij + kij into Eq. I(41|l to obtain 
5h(l ab ): 



Sh(l ab ) = iV7 mat — f k a - —kiJbakbi 
m \ Zir 



x ( kj k - —kjJbahk j E'^ (51) 



TVS 2 , 

= const -7 m at- E tj 

zirm J 



(52) 



^■(kia^bakbkkjk ~\~ kj a lb a kbkkik) 



The const= N^ ma ,t^E' ij 5ikkukjk does not depend on l ab . 
A term of order £^ -(Z a j,) 2 has been neglected. 



^dc . ( kd a l ca -\- k ca lda 

Air V 



(59) 



This can now be used to express the elastic energy in 
terms of l a b'- 



C4last — nKjkl-7— (ki a lja + kj a li a ) — (klblkb + kkbhb)V . 
Z 47T 47T 



(60) 

Some manipulation of indices using \abcd — ^-bacd — ^cdab 
leads to the more simple form 



V 



(61) 



C. Minimizing the total energy per unit cell 



B. Elastic energy 



The elastic energy is given as [38 



C4last — -^KjklCijeklV , 



(53) 



where Xijki is the elastic modulus, V the volume consid- 
ered, and 



dx-j 



dxi 



(54) 



is the strain tensor, where Ui is the displacement of a 
volume element at some location xi. For i = j, Uij repre- 
sents the relative change of length in a^-direction, and for 
i 7^ j, it represents the change of the right angle between 
lines originally pointing in Xi and Xj direction. 

To express the elastic energy in terms of lij, we note 
that the location x b of a point of the direct lattice can be 
expressed as a linear combination of the primitive lattice 
vectors 



(55) 



Summing up the contributions, we find for the energy 
change per unit cell (leaving out the constant terms) 



1 



Nh 2 



U = Idet^i^l-^-^XijklkiJjaklblkb - 7mat^— (62) 

^E i j[ki a li )a ki ) kkjk + kj a lb a 

kbkhk) 

were V — |det(2y)|, the volume of a unit cell, and 
N = N e<u the number of valence electrons per unit cell 
have been inserted (not to be confused with N eiTn the cor- 
responding number per molecule). The inner-shell elec- 
trons are assumed not to influence the crystal geometry. 
A minimum is found, when 



dU 

Oimn 



= 0. 



(63) 



After some manipulation, the derivative can be expressed 
as 



dU 



|detfa)| - N e , u h 2 

T~n ^imklKinKlblkb ~ 7mat — (Oi) 

47r z Zirm 
^•E^Aki n k m kkjk 4- kj n k m kkik) = . 
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We denote Et^s — |(J3y + E'^) the symmetric part of 
the tensor E'^ . The last equation can be simplified a bit 
by multiplying with l pn : 

^h 2 N etU 

{pj) mk ik= (65) 



\mklklbhb ~ 7n 



For solving this for Z, we need the inverse fj. a bki (called 
the compliance tensor) of X abcd , defined by 

Habkl^abcd = Skc$dl ■ (66) 

Multiplying Eq. (|65(l with fidepm gives 
47rJV P „ ft 2 



kebldb = 7n 



EfpjjUdempkmkkjk (67) 



|detftj)| 

and a further multiplication by l es yields 
2h 2 N e u , 

Ids = 7mat^|^^^ -j| E^pj^/Idempleskmkkjk ■ (68) 

The strain tensor can be calculated from this result using 
Eq. (El as 



e dc — BdcpjE'( p j) 



(69) 



with 



2N e u h 2 

Bdcpj — Tmat T~i TTr \T f-^dcmpkjnkkjk • C^) 

m\aet{lij)\ 

This has been simplified by using Eq. (|46|) . Since for Eq. 
1)69(1 . this is multiplied with the symmetric E [j p y only 

the part S^cpj = \{Bdcpj + Bdcjp) that is symmetric in j 
and p 

Bdcjp — '"/mat ~7\ 77 (fJ>dcmpkmkkjk ~t~ l^dcmjkmkkpk) 



' m\det{lij)\ 
is used. The resulting strain tensor 

*-dc $dcpj Epj 



(71) 
(72) 



is given by the 3x3 Lorentz violation tensor ■ and a 
tensor Bdcpj , which gives the sensitivity of the material 
geometry change to Lorentz violation in the electron sec- 
tor of the SME. This is the desired result of this section. 



1. General properties of the sensitivity tensor Bdcpj and 
conventions 



The sensitivity tensor has the symmetries 

Bdcjp — B c( jjp — $dcpj ■ (^^) 

In general, however, Bdcjp 7^ Bjpdc- It has, therefore, 
at most 36 independent elements, the number of which 



is reduced for a symmetric crystal. For a compact pre- 
sentation of the material specific results in the following 
sections, we will arrange these into a 6 x 6 matrix, that 
allows to express Eq. (|69l) as a 6 dimensional matrix 
equation 

e r = BteE' s . (74) 

We therefore arrange the six independent elements of edc 
and as the vectors 

— {pxx : &yy s Czz j ^-yz > &zx : &xy ) > (^^) 
E T = { E 'xx> E yy> E 'zz> E yz> E zx> E 'xy) ( 76 ) 

[the capital greek indices run from 1 ... 6] and define the 
sensitivity matrix 



B 



(77) 

The factors of 2 account for the double-counting of the 
non-diagonal elements of E '( p j-* in the tensor equation Eq. 



2. Sensitivity matrix for isotropic materials 

Let us first consider isotropic materials that have no 
preferred crystal orientation, i.e. crystals of cubic struc- 
ture and non-crystalline (fused) materials which consist 
of a large number of small crystals oriented statistically. 
Cubic materials have one single lattice constant a; the 
matrix of the primitive direct lattice vectors is given by 
Uj = aSij. According to Eq. I|46|l . the matrix of the re- 
ciprocal lattice vectors is given by fey = (2Tr/a)Sij. In 
the appendix, it is described how to obtain the compli- 
ance constants [i a bcd from the elasticity constants that 
are tabulated for various materials in the literature, e.g. 
[39l ]. Inserting into Eqs. (|71I77|I . we obtain the sensitivity 
matrix B. It is of the structure 



/Biui 


01122 


01133 


201123 


2B1131 


201112 




Briw 


02222 


02233 


202223 


202231 


202212 




Bzzn 


03322 


03333 


203323 


203331 


203312 




B23 11 


02322 


02333 


202323 


202331 


202312 




£>3111 


03122 


03133 


203123 


203131 


203112 




V B1211 


01222 


01233 


201223 


201231 


201212 


J 



(Bn 


012 


B12 













Bu 


B n 


B12 













B12 


Bl2 


B12 






















Bu 






















Bu 







V 














Bu 


) 



B 



For cubic crystals, the non-zero values 



(78) 



Bu 
B12 
Bu 



Z/(a 2 )Cii , 
C/(a 2 )Ci2 , 
Z/{2a 2 )Cu 
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where 



£ = 7mat 



8n 2 N e , u H 2 
m\det(kj)\ 



(79) 



and Crs are the elements of the 6x6 compliance ma- 
trix, Eq. (|lj9fl . From symmetry arguments, this is also 
the structure of the B matrix for non crystalline materi- 
als without a preferred orientation. The elements of this 
matrix for some cubic and/or fused materials are given 
in table G] The values for fused quartz and sapphire have 
been calculated from the values of the crystalline ma- 
terials (to be calculated below) as averages over crystal 
orientations. 



TABLE I: Elements of the sensitivity matrix for fused and/or 
cubic materials, fq denotes fused quartz, fs fused sapphire, 
C denotes diamond. Materials for which three elements of B 
are given are isotropic; the coefficients should be inserted into 
Eq. 1781 . The other materials are trigonal; the coefficients 
for these are to be inserted into Eq. 18211 . 



Mat. 


7mat 


Bii 


Bl2 




8m 




B 33 


B41 


B44 


fq 


0.38 


0.77 


-0.09 












0.57 


fs 


0.29 


0.06 


-0.01 












0.05 


Si 


0.50 


2.51 


-0.70 












2.05 


c 


1.16 


5.77 


-0.59 












5.35 


A1 2 3 
Si0 2 


0.29 
0.38 


0.14 
1.41 


-0.01 
-0.07 


-0.00 
-0.06 


0.02 
0.35 


-0.02 
-0.14 


0.01 
0.44 


0.01 
0.25 


0.04 
0.41 



3. Sensitivity matrix for trigonal crystals 

Quartz and Sapphire are of trigonal structure and are 
frequently used in cavity experiments. Therefore, we also 
consider the trigonal case here. The matrix of the prim- 
itive direct lattice vectors can be chosen as 



hj — 



( a/2 a/2 
V3a/2 -V3a/2 
^ c 



(80) 



where a and c are the two lattice constants. We calculate 
the inverse kij\ the product 



k fcjk 




(81) 



turns out to be a diagonal matrix. Trigonal crystals 
have six independent compliance constants and two lat- 
tice constants, which makes 8 independent components 
for the B-matrix. It has the structure 



B = 



( fill 


B12 


B 13 


B U 










3£ 12 




B13 


-Bl4 










B31 


§B 31 


B33 













B41 







£44 






















B55 


3B41 




V 











B55 


B 66 


) 



(82) 



Bl4 

B31 
B 33 
B41 
B 44 



£[l/(3a 2 ) + l/(2c 2 )]<7 14 . 
2£/(a 2 )C 13 , 

e/( c 2 )c 33 , 

£/(a 2 )C 14 , 

£[l/(6a 2 ) + l/(4c 2 )]C 4 4 ■ 



Crs are the elements of the compliance matrix given in 
Eq. I|B9(I . Numerical values of Brs for quartz and sap- 
phire are given in Tab. [I] The matrix B is not sym- 
metrical; the elements of the first column are generally 
the highest in this matrix, i.e., the geometry change of 
trigonal materials is most sensitive to the xx element of 
the Lorentz violation parameters E'ujy This is because 
the direct lattice vector components in ^-direction are a 
factor y3 smaller than the y components. Hence, the 
wave function of the electrons oscillate faster in x direc- 
tion, i.e. the (p x ) momentum component is larger. Since 
the influence of Lorentz violation is given by the (piPj) 
matrix element, this means a higher influence of the x- 
component of the Lorentz violation coefficients E',,-s. 

High elastic constants decrease the values of B so that 
crystals of high stiffness (such as sapphire) should show 
lower values of B. However, in some cases (particularly 
diamond) , this is outweighed by small dimensions of the 
unit cell (that imply high momentum expectation values 
due to short period of the electron wave functions) and 
a large number of electrons per unit cell. 



with 



EXPERIMENTS 



B55 = Ba + 2^41 , 
Bqq — 2^11 ~ ^12 • 
The matrix elements are explicitly 

B u = 2£/(a 2 )C n , 
B12 = 2£/(3a 2 )(7 12 , 
B13 = £/(c 2 )C 13 , 



Here, we discuss the application of our results to ex- 
tract limits on Lorentz violation in the electrons' equa- 
tion of motion from experiments. 



A. Lorentz violation signal in cavity experiments 

As discussed in the introduction, Lorentz violation may 
affect the resonance frequency of a cavity v through a 
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variation of c and L: 



5c 

CO 



SL 



(83) 



Here, vq,Cq, and Lq are quantities in absence of Lorentz 
violation. For cavity experiments, we have to consider 
both Lorentz violation in electrodynamics as well as 
in the electrons' equation of motion. The influence of 
Lorentz violation in electrodynamics leads to a frequency 
change Semv that is mainly caused by a variation of the 
velocity of light, <5emc 0, 0, and a small material- 
dependent contribution due a length change S e -L, that 
is usually negligible [l9j| . Lorentz violation in the elec- 
trons' equation of motion affects solely L. This leads to 
a frequency change 5 e -V, so that 



8v = Semv + 5 e -i> . 



(84) 



For example, if the cavity axis is parallel to the z-axis of 
the crystal, 



8 e -i> SL Z 



e3 • 



(85) 



From Eq. O, e 3 = B 3S E' S . 

B. Limits from previous experiments 

Cavity experiments have been performed repeatedl y as 
tests of Lorentz violation in electrodynamics 0, TToL TTlL 
H3 . ITsl IbH . ITfij . It is interesting to estimate the level 
of the limits on c M „ resulting from these experiments. 
Since none of them were done with a setup optimized for 
obtaining separate bonds on Lorentz violation in electro- 
dynamics and in the electrons' equation of motion (see 
below), we have to use an assumption that simplifies the 
analysis, so we can work with the published data only. 
To obtain a sharp upper limit, we have to compare two 
experiments of high precision that used different cavity 
materials. 

Muller el al. [l^ performed a Michelson-Morley exper- 
iment using two sapphire (s) cavities subject to Earth's 
rotation. The cavity axes were parallel to the crystals 
c-axis. Since for such cavities, the elements of Bg S giv- 
ing the sensitivity to Lorentz violation in the electronic 
sector are relatively low, we neglect the effect on Lorentz 
violation in the electrons' equation of motion for this ex- 
periment, i.e., we view the experiment as a pure test 
of Lorentz violation in electrodynamics. Due to Earth's 
rotation with the angular frequency ~ 27r/23h56min 
and Earth's orbit with fi© = 2ir/l year, Lorentz violation 
in electrodynamics leads to a time-dependency of the fre- 
quency difference between the cavities that has Fourier 
components at 6 frequencies that are linear combinations 
of and O ffi . The experiment gives individual bounds 
on the amplitudes of these Fourier components. At 2u>^ 
and 2(oj0 ± fi®), the bounds are below 4 x 10~ 15 . 



Brillet and Hall Q used a single fused quartz (fq) 
cavity (actually, "ultra- low expansion" glass ceramics, 
ULE) on a turntable rotating at uj t . Frequency measure- 
ment was accomplished by comparison to a stationary 
CH4 standard. Since in a single-cavity setup, the same 
Lorentz violation in electrodynamics leads to half the fre- 
quency variation compared to a two-cavity setup, the ex- 
periment of Muller et al. excludes signals from Lorentz 
violation in electrodynamics larger than 2 x 1CP 15 for 
this experiment. We may thus view this experiment as a 
measurement of 



3H 



(86) 



The indices S denote components in the cavity frame of 
reference. In accordance with [TtI fljl |23| , we define the 
x axis of the laboratory frame as the north-south axis, 
the y axis as the east- west axis, and the z axis as point- 
ing upwards. The turntable rotated in the ccy-plane 0. 
We define the cavity frame z axis parallel to the z axis, 
the x axis parallel to the cavity axis. A calculation of 
the hypothetical signal starts from transforming c M „ as 
given in the sun-centered standard frame into the lab- 
oratory frame, as described in 0, [LsL |23| . A further 
rotation around the laboratory z axis gives the quanti- 
ties Cpp in the cavity frame, from which Ej-. follows from 
Eq. JSJ. Due to the rotations, the hypothetical signal 
becomes time-dependent and is given by (assuming that 
at T ffi = 0, the cavity axis coincides with the X axis, and 
neglecting terms proportional to Earth's orbital velocity 

Ad ~ 10- 4 ) 
Sv 

— = C"(2,0,0)cos2w t T e + S'(2,2,0)sin(2a; t + 2a;e)T e 
v 

+C(2, 2, 0) cos(2^ t + 2 Wffi )T ffi + A + 0{p 2 m ) , (87) 

where A denotes Fourier components at other frequen- 
cies, for which no experimental results are published in 
0. The sine component at 2u> t vanishes. The coefficients 



C(2,0,0) 
5(2,2,0) 

C(2,2,0) 



1 



- B{!) sin' X b(cxx + c yy - 2c zz ) 

-0.13 (c X x + c Y Y - 2c zz ) , 
C0S ^ iB h_ B f i)c(XY) 

0.58 C(xy) , 

l - CO s^{B[l-B[«){c YY -cxx) 

0.29 (cyy-Cxx). 



B fq 



Xb denotes the colatitude of Boulder, \b — 50° 
and #(2 are given in Tab. [I] From the experiment, an 
amplitude at 2u>t of ~ 2 x 10 -13 is found, attributed to a 
slight tilt in the horizontal alignment of the turntable. In 
principle, a hypothetical Lorentz violation signal at 2uj t 
cannot be separated in the analysis from a signal caused 
by such tilt. However, if we consider an exact cancellation 
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between a strong Lorentz violation signal and a strong tilt 
signal improbable, we obtain a rough order-of-magnitudc 
bound on \cxx + cyy — 2czz\ at a level of about 1CU 12 . 
The measured upper limit on the amplitude of the com- 
ponent of the hypothetical signal at 2(u> t + w^) (which 
can be distinguished from a tilt generated signal in the 
Fourier analysis) is 4 x 10 -15 . Adding in quadrature the 
maximum possible contribution from Lorentz violation in 
electrodynamics according to the result of |l5j, 2 x 10 -15 , 
we obtain ^5(2, 2, 0) 2 + C(2, 2, 0) 2 < 4.5 x 10" 15 . This 
gives the limits 

\c { xy)\ < 8x 1(T 15 , (88) 
\cxx - c YY \ ^ 1.6 x 1(T 14 . 

Given the small magnitude of these values would in- 
deed lead to negligible contributions to the experiment of 
Miillcr et al, so our above assumptions seem reasonable. 

C. Optimized setups 

To obtain clean separate bounds on Lorentz violation 
in electrodynamics and the electrons' equation of motion, 
without using the above assumptions, a dedicated exper- 
iment is desirable. One could, for example, compare the 
resonance frequencies Vf q and v s of a cavity made from 
fused quartz and one from crystalline sapphire, with the 
cavity axis parallel to the crystal z axis, for example. 
This seems to be a realistic scenario, since such cavities 
have been used in exp eriments and proved to be of high 
stability H [H [H El E M, E3 The signal for Lorentz 
violation would be the frequency difference v s — Vf q , If 
a parallel arrangement of cavities would be chosen, the 
influence of Lorentz violation in electrodynamics could 
be eliminated, and a clean bound on some components 
of E'- k could be extracted, 

^l^L = m(i - B S ^)E'^ . (89) 

Here, v ~ v s ~ v q is the average frequency. An accu- 
racy level of below one part in 10 15 in frequency com- 
parisons of cavities made from quartz |40j and sapphire 
[15| has been demonstrated in the laboratory. Thus, 
placing bounds of a few parts in 10 15 on the compo- 
nents of E' JK that dominate the signal seems feasible. A 
contribution of the time-components cmj) to the signal 
arises if one takes into account the laboratory velocity 
given by the velocity 0® ~ 1CP 4 of Earth's orbit and 
< /3l ^ 1.5 x 10~ 6 due to Earths rotation (depending 
on the geographical latitude). The Lorentz transforma- 
tions between the sun-centered inertial reference frame 
and the laboratory frame lead to additional Fourier com- 
ponents of the signal that are proportional to C(oj) and 
either /3® or (3l- In a Fourier analysis of a sufficiently 
long timetrace, the Fourier components can be resolved 
and individual limits on almost all components of c^ v 
(only cqo does not lead to time-dependent signals to first 



order in /3 e or fti) can be expected, at or below about a 
part in 10 15 for the dominating parameters and to about 
a part in 10 11 for the parameters that are suppressed by 
/3®. Future space experiments, for which a resolution of 
the frequency measurement of up to 10" 18 is projected 
HHHI). might bound the dominating components of Ej k 
at the 10 -18 level, and the suppressed components of c^ v 
at the KT 14 level. 

Instead of using different cavity materials, two cavi- 
ties made from the same crystalline material, but having 
different orientations of the cavities with respect to the 
crystal axes might be used. Using, e.g., quartz, a compar- 
ison between a cavity fabricated such that the cavity axis 
is parallel to the crystals c axis to one having its cavity 
axis perpendicular to the c axis would provide a relatively 
high sensitivity given by B\~ and B\~ . Another possibil- 
ity would be two orthogonal cavities within a single block 
of crystalline material. Such an arrangement might be 
favourable for eliminating parasitic effects, like thermal 
expansion and vibration. However, in such an experiment 
the frequency change <5em^ due to Lorentz violation in 
electrodynamics would not drop out, thus complicating 
the analysis. Another configuration would be the com- 
parison of a cavity against an atomic clock, such as a 
cesium clock and/or a hydrogen maser. This scenario is 
interesting, since it is projected for the OPTIS satellite 
|2lj . This sattelite is projected to carry three cavities 
orthogonal to each other, so it is probably possible to 
separate the electrodynamic terms from the electronic 
ones in an analysis of the complete data, that consists of 
all frequency differences between the three cavities and 
the atomic clock(s). Such an analysis would also have 
to take into account a possible Lorentz-violating shift of 
the atomic clock frequencies, which could give the exper- 
iment sensitivity to additional parameters. 



VI. SUMMARY AND OUTLOOK 

We have calculated the change of the geometry of crys- 
tals that is caused by Lorentz invariance violation in 
the fermionic sector of the extended standard model. 
The length change is caused by a modified kinetic en- 
ergy term (djk + 2E'- k ){l/2m)pjPk that enters the ki- 
netic energy term of the Hamiltonain for the free elec- 
tron. E'j k = —Cjk — ^cooSjk, where c^ v is a Lorentz ten- 
sor originating from the standard model extension. The 
calculation proceeds using a Bloch ansatz for the wave 
function of the valence electrons with the lattice periodic 
function given by a Fourier series. The crystal adjusts its 
geometry such as to minimize its total energy. In that 
way, Lorentz violation in the electrons' equation of mo- 
tion affects the length of an electromagnetic cavity that is 
made from the crystal, and thus the resonance frequency 
of a cavity made from the material. As a main result of 
this paper, there is thus a method to measure the in 
cavity tests of Lorentz violation. 

Comparing cavities made from different materials, it 
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is possible to separate the effect connected to c M „ from 
Lorentz violation in electrodynamics, that also affects the 
resonance frequency of cavities. Under some assumptions 
that help to separate the electrodynamic and the elec- 
tronic terms, already performed experiments 0, IbH Il5j 
imply constraints on c^xy) an d cxx — cyy at the 10~ 14 
level, and on cxx + cyy — 2czz at the 10 -12 level. This 
is to our knowledge the only present experimental con- 
straint on the components of c^ v . We discuss possible 
setups for experiments that can obtain separate bounds 
without using these assumptions, and obtain results on 
more components of c^ v . Future experiments on Earth 
and in space promise increased sensitivity up to a part 
in 10 18 . 

In the appendix, we briefly discuss the case of 
spin-polarized matter. An additional contribution to 
the length change arises from a spin-dependent term 
Fjki^PjPkP 1 also originating from the standard model 
extension. This allows to deduce limits on Fj fc; from ex- 
periments using a spin-polarized cavity material, at least 
in principle. 

Our model of the solid state could be improved by us- 
ing material specific values for the Fourier coefficients 
of the single electron wave function. Since our model 
Fourier coefficients already satisfy the symmetry require- 
ments for a realistic wave function, this might result in 
relatively minor corrections for the length change. A 
check (maybe for a simple material), however, might be 
worthwhile. Most importantly, however, a dedicated ex- 
periment will be performed to obtain more complete and 
/ or stronger limits on Lorentz violation in the electrons' 
equation of motion. Cavities made from crystalline sap- 
phire and fused quartz are ready to be implemented. 



1. Hamiltonian 

If the average of the spin expectation values is non- 
zero, the spin-dependent terms contribute to the Lorentz- 
violating correction to the hamiltonian. We still assume 
a vanishing average helicity, Eq. I|22[l . for all electrons. 
Thus, 

JV -, N 

(Sh) = mc i B' :j J2((a a y) + -Ei :i J2((Pa)i(Pa) j ) 

a—1 a—1 

1 N 

+- F kiT,((p^M^y) ■ ( ai ) 

a=l 

We assume that a fraction rf of the total N electrons have 
their spin \ polarized parallel to x l . The other electrons 
are assumed to be unpolarized. The average of the spin 
expectation values is then given by 

a 

For the last term of Eq. (|A1|) , 

N N 

E<(Pa)i(P-)*(^)') = 2^E«^M^> ■ ( A3 ) 
a—1 a—1 

Therefore, 

(Sh) = mc 2 ^-B'-rf (A4) 
2 J 

1 / 1 \ N 

+ (E[, J +-Fl 3kV k )Y J ({PaUPa) ) ■ 
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APPENDIX A: SPIN-POLARIZED MATERIALS 

For a spin-polarized material, an additional contribu- 
tion to the geometry change arises from a spin-dependent 
term of the non relativistic single-electron hamiltonian of 
the SME, Eq. ©, that is given by Fj^i ■ That means, 
from experiments using cavities made from spin-polarized 
materials, a limit on Fj kl can, at least in principle, be de- 
duced. This is interesting, since many degrees of freedom 
of F'j kl are not yet fixed experimentally. In this appendix, 
we estimate the effect and the level of sensitivity that can 
be expected for such an experiment. 



2. Geometry change 

The term mc 2 N B'^rf /2 contained in the hamiltonian is 
independent from the crystal geometry and does, there- 
fore, not lead to a geometry change. The second term 
that is proportional to the average of (piPj) over all 
electrons, however, leads to a geometry change, that 
can be calculated in analogy to the discussion in the 
main parts of this paper. We can overtake the result, 
Eq. I|tj9fl . for the geometry change if we replace E'^ by 

E 'ij = E ij + iKjk^- 

The sensitivity of the cavity geometry to F'- kl is thus 

given by Brs as well as rf. The magnitude of the lat- 
ter can be estimated as the ratio of the number of spin- 
polarized electrons ns, u to the total number of electrons 
N e>u per unit cell, \q\ = riB,u/n eyU . In a saturated ferro- 
magnetic material, e.g. iron at a magnetic field of 1.7 T, 
Ub,u ~ 2.2 spins are polarized per unit cell, so \r}\ is 
of order unity. (Note that ub,u can be as high as ~ 10 
for Dysprosium.) Therefore, the sensitivity of the cav- 
ity length (and thus, its resonance frequency) to Fj kl is 
comparable to the sensitivity to E'- k . 



14 



3. Possible experiments 



and 



If the cavity is made from a spin-polarized solid, i.e., 
a magnetized ferromagnetic material, the cavity length 
would depend on E'^ + F-j k r) /2. That means, from a 
measurement of the resonance frequency of such a cavity, 
a limit on F[- k could be derived, provided that separate 
limits on E[- are known from previous experiments us- 
ing one of the methods discussed above. However, note 
that the systematics of such an experiment are largely 
unknown. The selection of materials suitable for build- 
ing stable cavities is a highly nontrivial discussion of the 
experimental systematics, some of which are far from ob- 
vious. Because of the manifold effects connected to mag- 
netization (e.g., magnetostriction), an experiment using 
a magnetized cavity could suffer from various systematic 
effects, so our discussion is a bit speculative. A theoreti- 
cal complication is that practical ferromagnetic materials 
are usually alloys (such as AINiCo), whereas the theory 
presented above is directly applicable for crystals only. 

For such a measurement of F^ k , one could use a cavity 
made from a permanent magnetic material. The direc- 
tion of the spin polarization with respect to the cavity 
axis determines the components of F^ k which dominate 
the experiment. A rotation of the cavity would modu- 
late the .F/^-induced frequency shift. The corresponding 
time dependency of the cavity resonance frequency would 
be the signal for a non-zero F!- k , It could be beneficial to 
use solely the Earth's rotation to avoid possible system- 
atics associated with a magnetized cavity rotating in the 
Earth's magnetic field. Magnetic shielding will probably 
also be necessary. If the frequency stability of a cavity 
made from a suitable magnetized material would be of 
the same order as the stability achieved with quartz or 
sapphire cavities, limits on F^ k of order 10 -15 could be 
achieved. 
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This is called the Voigt convention in the literature |37j . 
pp. 604-609. The matrix S is symmetric, it thus contains 
at most 21 independent elements. The symmetry of the 
crystal reduces the number of independent elements. For 
example, the matrix for cubic symmetry has three inde- 
pendent elements: 



/ Six S12 S12 
S12 an S12 
S12 S12 Sn 
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) 544/ 
For trigonal symmetry, there are six: 



S = 





S12 


Sl3 


5*14 








\ 


S12 


Su 


Sl3 


— S14 










5*13 


Sl3 


S33 













5*14 


— S14 





S44 






















S44 


S14 




V 











Sl4 


"5*66 


/ 



(B7) 



where 



2(Sii — S12 



(B8) 



APPENDIX B: NOTATION CONVENTION IN 
ELASTICITY THEORY 

For obtaining the material specific values of the sensi- 
tivity tensor B a bcd, the compliance tensor n a bcd has to be 
known. The relation between stress Oij and strain e^c is 
given as 



(Bl) 



In engineering, it is common to replace this by a six- 
dimensional matrix equation 



a = S • e . 



where 34], p. 445 



a 
e 



iPxxi ®~yyi 0~zz: &yz: &zx: ®xy) 3 
(&xx: &yyi &zzi ^^yz: *2&zx> 2e X y) 1 



(B2) 



(B3) 
(B4) 



By inverting S one obtains the compliance matrix C that 
satisfies S • C = 1, where 1 is the six dimensional unit 
matrix. C enters the equation e = C • a between the 
stress and strain 6- vectors Eqs. l|B3IB4|l . On the other 
hand, the compliance tensor /j, a bcd enters the relationship 
between the stress and strain tensors, e a b = fJ-abcdO'cd- 
Therefore, the compliance tensor is related to the com- 
pliance matrix by 



C = 



/ /^llll A*1122 M1133 

^2211 A*2222 ("2233 

^3311 ^3322 M3333 

2/13211 2/13222 2/i3233 

2/i3111 2/i3i22 2/13133 

V 2/11211 2/11222 2/11233 



The compliance and elasticity 
tals have the same symmetry; 
symmetry of the compliance 
one of the elasticity matrix, 



2/11123 2/1H31 2/11H2 \ 

2/12223 2/12231 2/12212 

2/^3323 2/I3331 2/I3312 

4/13223 4/13231 4/13212 

4/13123 4/l3i3i 4/13H2 

4/11223 4/11231 4/11212 / 

(B9) 

matrices for cubical crys- 
for trigonal crystals, the 
matrix is similar to the 
with the exception that 
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Ce6 — 5 (Cn — C12). The compliance tensor elements 
fi a bcd a re obtained from the tabulated elements of the 
elasticity matrix S by inverting the elasticity matrix and 
reading of the tensor elements from Eq. I|B9|1 . 



where \ is the geographical colatitude, and ~ 
27r/23h56min Earth's rotation angular frequency, leads 
to the tensor c M „ as expressed within the laboratory 
frame. Another rotation around the z axis using the 
rotation matrix 



APPENDIX C: SIGNAL COMPONENTS FOR 
LABORATORY EXPERIMENTS WITH 
TURNTABLE 



cos ujttt sin cotU 
Rt — I — sin U t t t cos ujttt 
1 



(C3) 



Here, we give the full signal components caused by 
Lorentz violation in the electrons' equation of motion 
in the laboratory frame for a cavity rotated, using a 
turntable, at an angular frequency u> t , assuming a ma- 
terial of trigonal or higher crystal symmetry, i.e., the 
sensitivity matrix B is of the form Eq. (|82|) or simpler. 
The rotation axis is fixed to point vertically. We use a 
turntable time scale U defined such that tt = at any 
one instant when the cavity is pointing in the x direction 
of the laboratory frame. 

We use two reference frames, one sun-centered celestial 
equatorial reference frame and one laboratory frame. As 
defined in [2^, the sun-centered frame has the X axis 
pointing towards the vernal equinox (spring point) at Oh 
right ascension and 0° declination, the Z axis pointing to- 
wards the celestial north pole (90° declination) and the 
Y axis such as to complete the right handed orthogonal 
dreibein. Earth's equatorial plane lies in the X — Y plane; 
its orbital plane is tilted by 77 ~ 23° with respect to the 
latter. The time scale T = when the sun passes the 
spring point, e.g., on march 20, 2001 at 13:31 UT. 
The laboratory frame has the x axis pointing south, the 
y axis east, and the z axis vertically upwards. The lab- 
oratory time scale T$ = when the y and the Y axis 
coincide. 

The signal derivation starts from the symmetrized ten- 
sor c^ v ) given in the sun-centered celestial equatorial ref- 
erence frame, which is suitable for expressing the tensor 
because it is inertial on all time-scales involved in ter- 
restrial experiments. To C( Mt/ ) , we first apply a Lorentz 
boost to first order in /3® ~ 10 -4 , the velocity of Earth's 
orbit 



' sin T 

— cost? cos f2 ffi T 
y — sin r\ cos ^1®T 



(CI) 



where Q,® ~ 27r/lyr is the angular frequency of Earth's 
orbit. We neglect the smaller velocity 0</3l<1.5x 
10~ 6 due to Earth's rotation in order not to complicate 
the signal components below further. Subsequently, ap- 
plication of the rotation matrix 

cos x cos u)($T(£ cosx sin w®T® — sinx 
R = I — sin a;®!® cosw®T0 
sinx cos w e T e sinx sin w®T e cosx 

(C2) 



leads to the quantities within the rotating turntable 
frame, which are then decomposed according to Eq. JHJ. 
The time scale t t = when the cavity axis is parallel to 
the x axis. Insertion of the results into Eq. I|74|l gives 
the cavity length change, and thus the frequency change 
that is given below. 

For compact notation, we define the abbreviations 



u)(a, 6, c) = auj t + 6u>© + cQq 
4>(a, b, c) = auj t t t + 6w e T e 4 



(C4) 

cft ffi T. (C5) 



We give the signal components for the two most inter- 
esting cavity constructions: A cavity with the resonator 
axis pointing parallel to the crystals 's c axis (which is cur- 
rently the most familiar cavity type), and a cavity with 
the resonator axis parallel to the crystal's a or b axis 
(which gives sensitivity to the C( i) components to first 
order in f3<$, the Earth's orbital velocity). The signals 
are expressed as a Fourier series 

— = C(0,0,0) + ^[5(a,6,c)si n(? !)(a,6,c) (C6) 

a.b.c 

+C(a, b, c) cos 4>{a, b, c)] 

with coefficients S(a,b,c) and C(a, 6,c); the dc compo- 
nent C(0, 0, 0) is not included in the equations below, as 
it is not measurable. 



1. Cavity axis parallel to a axis 

We use the abbreviations Ba = Bn + B12 — 2/3i 3 and 
Bb = Bu — B\2- The signal consists of 18 frequencies 
w(a, b, c) with 

C(0,1,0) = -B A c (X z) cosx sinx, 

5(0,1,0) = -B A C(YZ) cosx sinx, 

C(0,2,0) = -(l/4)B A (cxx -cyy) sin 2 x, 

5(0,2,0) = -(l/2)/3 A c (xy) sin 2 X , 

C(l,-2,0) = 2/3 14 c (xy) cos|sin 3 |, 
5(1,-2,0) - /3 14 ( Cxx - Cyy ) cos ! sin 3 |, 



C(l,-1,-1) = -/3 14 /3esin 2 |sin|( 

• V\ 
+C(TZ)Sm-j , 



C(TY) COS - 



16 



5(1,-1,-1 
(7(1,-1,0 
5(1,-1,0 
C(l,-l,l 



5(1,-1,1 

(7(1,0,-1 
5(1,0,-1 
(7(1,0,0 
5(1,0,0 
(7(1,0,1 
5(1,0,1 

(7(1,1,-1 



5(1,1,-1 
(7(1,1,0 
5(1,1,0 
(7(1,1,1 

5(1,1,1 
(7(1,2,0 
5(1,2,0 

(7(2,-2,0 
5(2,-2,0 
(7(2,-1,0 
5(2,-1,0 
(7(2,0,0 



-2 B i4Pffi c (TX)Sin -sm?7, 
BuC(yz) (1 + 2 cos x) sin 2 - , 
B 14 c {xz) {l + 2cos X ) si n 2 |, 



2?i 4 /3 e cos | sin 2 | 



: (TZ) 1 



-c (TF) sin - 

- 2 S i4/?e c (TX) sin 2 - sin i) , 

-(l/2)Si 4 /3eC(TX) cos ?7 sinx, 

(l/2)S 14 /3 e c (Ty) sinx, 

0, 

- (1/4)B U (c XX + cyy - 2c zz ) sin 2 X . 
-(l/2)0i 4 /3 ffi c (TX) cos sin x, 
-(l/2)2?i 4 /? e c (Ty) sinx, 

-Bi 4 /3e co s 2 cos 2 ( C ( TZ ) cos 2 
. »7 

-C(TF) sm - 

1 k a 2 X . 

-2»i4Pe c (TX)Cos -sin 77, 

BuC(yz) cos 2 I (2 cos x - 1) , 

-BuC(xz) cos 2 I (2 cos x - 1) , 

n n 2 X • V ( V 
Ci 4 /?e cos - sm - [C(ty) cos - 

+C(tz) sm - 

1 k a 2 X . 

-2#i4PeC(TX)Cos 2 sm7 '' 

(l/2)i3i 4 C(xy)(l + cosx) sinx, 

r ( \ 3 X ■ X 

-B u {c X x - cyy) cos ^-s 111 ^ 

^B B (c X x - cyy) sin 4 |- , 



K> • 4 X 

-B B C(XY) sm 



2 ' 
X • 3 X 



26 B C(xz) cos sm" ^ 
-2B B C(y Z) cos ^ sin 3 - , 
-(l/4)B B {c xx + c YY - 2c zz ) sin 2 X ■ 



5(2,0,0) 
(7(2,1,0) 

5(2,1,0) 

(7(2,2,0) 
5(2,2,0) 



= -26 



3 X 



bC(xz) cos - sm - , 
-2B b C(yz) cos 3 ^ sin ^ , 
^B B (c X x - cyy) cos 4 ^ , 



= £ 



4 X 



scpcr) cos - 



The signal for a cavity parallel to the crystals b axis can 
be obtained from these equations, if the x and y axis are 
interchanged. 



2. Signal for a cavity parallel to the c axis 

The cavity is oriented with its axis parallel to the Xt 
axis. We introduce the abbreviations Be = — 2Z?3i + 
2?32+2?33 and Bd = #32—033- (For the trigonal case, Eq. 
(|S2*|l . fi 32 = /?3i/3, for isotropic materials, #32 = #31-) 
We have seven signal frequencies with the amplitudes 



(7(0,1,0 
5(0,1,0 
(7(0,2,0 
5(0,2,0 

(7(2,-2,0 

5(2,-2,0 

(7(2,-1,0 

5(2,-1,0 

(7(2,0,0 
5(2,0,0 

(7(2,1,0 
5(2,1,0 
(7(2,2,0 
5(2,2,0 



-Bcc(xz) cos x sinx, 
-BcC(yz) cos x sinx, 
-{\/A)Bc{cxx - cyy) sin 2 x • 
-(1/2)B C C(XY) sin 2 x, 

-^B D (c X x - c Y y) sin 4 ^ , 

K> . 4 X 

»ZJC(xy) sm - , 

X -3 X 



-2i3_DC (xz) cos - sin ^ 
2££,c (yz) cos I sin 3 | , 
{l/A)B D (c X x + cyy - 2c^z) sin 2 x ■ 



= 0. 



2Bdc (xz) cos 3 - sin ^ 



= 2B 



3 X . X 



dC(yz) cos 2 sm 2 
-^B D (c X x - c Y y) cos 4 |- 
-6 D c (X y) cos 4 I . 
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